Collapse of a semiflexible polymer in poor solvent 
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We investigate the dynamics and the pathways of the coUapse of a single, semiflexible polymer in a 
poor solvent via 3-D Brownian Dynamics simulations. Earlier work indicates that the condensation 
of semiflexible polymers generically proceeds via a cascade through metastable racquet-shaped, long- 
lived intermediates towards the stable torus state. We investigate the rate of decay of uncoUapsed 
states, analyze the preferential pathways of condensation, and describe likelihood and lifespan of 
the different metastable states. The simulation are performed with a bead-stiff spring model with 
excluded volume interaction and exponentially decaying attractive potential. The semiflexible chain 
collapse is studied as functions of the three relevant length scales of the pheno menon, i.e. , the total 
chain length L, the persistence length Lp and the condensation length Lq = ^ kBTLp/uo, where uq 
is a measure of the attractive potential per unit length. Two dimensionless ratios, L/ Lp and Lq/ Lp, 
suffice to describe the decay rate of uncoUapsed states, which appears to scale as {L/ LpY^'^ [Lo / Lp). 
The condensation sequence is described in terms of the time series of the well separated energy levels 
associated with each metastable collapsed state. The collapsed states are described quantitatively 
through the spatial correlation of tangent vectors along the chain. We also compare the results 
obtained with a locally inextensible bead-rod chain and with a phantom bead-spring model. Finally, 
we show preliminary results on the effects of steady shear flow on the kinetics of collapse. 

PACS numbers: 87.15.He, 36.20.Ey, 87.15. Aa 



1. INTRODUCTION AND MOTIVATION 

The conformation of individual polymer chains de- 
pends on the solvent properties 0, M 13 ■ In good sol- 
vent the monomers prefer being surrounded by solvent 
molecules, effectively repelling each other. This effect 
leads to a swollen coil conformation for flexible poly- 
mers in good solvent. In poor solvent, conversely, the 
monomers try to exclude the solvent and effectively at- 
tract one another, and a flexible chain forms a compact 
globule of roughly spherical shape to minimize the con- 
tacts between monomers and solvent. The dynamics 
of the coil-globule transition in flexible chains are well 
known 0, IM S ^-nd involve the formation of a pearl 
necklace and the gradual diffusion of large pearls from 
the chain ends 0, Q . 

Many polymers, however, exhibit substantial bending 
stiffness, i.e., they are semiflexible. Such polymers can 
be described by the worm-like chain model; examples in- 
clude biopolymers (F-actin, DNA) as well as synthetic 
polymers (xanthan, PBO, PPTA, kevlar). Semiflexible 
polymers in good solvents form open, extended struc- 
tures, while in poor solvent the effective self-affinity of 
the chain contrast this tendency and leads to chain col- 
lapse. The collapsed state configurations and the path- 
ways to their formation are the result of the interplay 
between two opposing potentials: the bending potential 
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related to the chain stiffness and the attractive potential 
due to the poor solvency of the environment. 

A compact globule is energetically unfavorable for 
semiflexible polymers because it involves highly bent 
states. The collapsed ground state is instead a torus, 
which reduces the monomer-solvent contacts without 
cau sing excessive bending penalty. Theoretical analysis 
[lol Hlf . Monte Carlo simulations ^3 and experimental 
evidences 0, have detailed the stability, the features 
and the packing of the collapsed tori. The dynamics of 
the collapse of worm-like chains in poor solvent have been 
investigated more recently 0, 0, , strongly suggest- 
ing a possible generic pathway for the collapse of semiflex- 
ible polymers, featuring a series of long-lived, partially 
collapsed, racquet-shaped intermediate states, before the 
eventual collapse to a torus. Such intermediate states 
form an energetically driven cascade of increasingly com- 
pact conformations with sharp transitions between their 
energy levels. 

In the present work, we verify this hypothesis and fur- 
ther explore the effect of the relevant length scales on the 
phenomenon. The natural length scale for the polymer 
stiffness is its persistence length Lp, defined as the ratio of 
the bending stiffness k to the thermal energy fcsT, which 
represents the length along the chain for which tangent 
vectors remain correlated in good solvent. The balance 
between attractive forces and bending forces can be con- 
sidered in term of the so-called "condensation length" 
Lq = ^kTLp/uo, where uq is the value of the attrac- 
tive potential per unit length. We choose as time scale 
for the collapse event the rotational diffusion time of a 
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perfectly rigid rod with L = Lp, rdiff,L = CL^/(72knT) 



where C is the transverse friction coefficient |18| . There- 
fore we use T^m^Lp as unit of time and fc^T as unit of 
energy in presenting aU the results of the present work. 
We investigate the sequences and the likelihood of differ- 
ent collapsed structures formed at various solvent quality 
(described by uq) and chain stiffness. 



2. SIMULATION DETAILS 

The dynamics of the continuous wormlike chain are 
described by a inertialess Langevin equation: 
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Equation n balances the hydrodynamic force exerted on 
the chain by the surrounding fluid with the force due to 
the global potential acting on the chain U = C/bend + 
t^conn + C4oiv, and the Brownian force per unit length r), 
with correlations {■n{s,t)r]{s' ,t')) = 2kBTC6{t - t')S{s - 
s'). C/bend = kTLp dsdu{s)/ds is the bending poten- 
tial related to the chain stiffness, C/com is the connector 
potential, UsoW describes the monomer-monomer inter- 
action, and the Brownian term accounts for the rate of 
fluctuating momentum transfer due to the random colli- 
sions of small solvent particles with the chain. 

The potential Uso\v chosen to account for the effect of 
poor solvent is an exponential attractive potential with 
an excluded volume repulsive term: 
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where uq defines the depth of the potential well, i?attr 
is the range of the attractive forces and a is the ra- 
dius of the excluded volume region. With this potential, 
uo = RattrUo- Different choices for the shape of the po- 
tential should not change the stable and metastable con- 
figurations for the wormlike chain in poor solvent, while 
they may modify the speed and the dynamics of the col- 
lapse. Our choice of the exponential potential is related 
to the physics of the condensation phenomenon. DNA 
and polyelectrolytes collapse because of charge inversion 
and counterion induced attraction 0. The attractive 
forces are therefore due to highly screened electrostatic 
charges, which typically show an exponential decay, mim- 
icked by the potential selected in this work. 

The Brownian Dynamics code used in this work solves 
in dimensionless, discretized form the Langevin equation 
of motion of a linear wormlike chain of N beads with 
positions Ri , . . . , Rat , connected by iV — 1 connectors 
of equilibrium length a = L/{N — 1) with unit tangent 



vectors = (Ri+i — Ri)/a. For most of the simula- 
tions reported in this work, we use quadratic springs as 
connectors, with a force constant equal to 100 kT, which 
ensures a variation of the connectors' length within 30 %. 
This choice allows for longer time steps without interfer- 
ing with the spring relaxation time. We use a midstep 
algorithm to compute the bead trajectories gener- 
ated by the equation of motion 
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Here, = is a bead friction coefficient, the bending 
force is: 

F'r" = ^/a E ■ (4) 

the elastic force due to the spring is 

Ff = H[(|u,_i| - a)u,_i - (|u,| - a)u,] , (5) 

the potential force related to the effective interactions 
between the beads of the chain — repulsive within the 
hard core, attractive at short distance is 
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where ARy = Ri — Rj, and F""'' is a random force. 

Another relevant issue in the simulations is the accu- 
rate description of the high curvature of the collapsed 
states. The discretized computational model for the 
wormlike chain retains its validity when subsequent con- 
nectors are very nearly parallel: to satisfy this condition, 
a large number of connectors within the characteristic 
length scale of the collapsed structure is needed. The 
typical size of the torus and of the racquet heads are of 
the same order of magnitude of Lp. A good accuracy in 
the description of their curvature requires a large number 
of connectors per persistence length; for all the simula- 
tions presented here, we use 20 connectors per Lp. 

We also performed simulations with a bead-rod model 
[21I , which ensures local inextensibility of the worm- 
like chain. The computational cost for this model is how- 
ever significantly higher (~ 10-50 times) than that for a 
bead-spring chain model, mainly because a smaller time 
step is required for convergence. While performing most 
of the simulations using the bead-spring model, we com- 
pared the results for a small number of conditions with 
the bead-rod model, as discussed in Section In that 
Section we also compare some of the results with that ob- 
tained using a bead-spring phantom chain, i.e., without 
excluded volume interactions. 

For simplicity, we use isotropic drag in the simulations, 
and we neglect hydrodynamic interactions between seg- 
ments of the chain. Including such interactions would 
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FIG. 1: Typical configurations of an ensemble of collapsing 
chains at an intermediate time {L/Lp — 8, uo — 1.25kT/a). 

likely give minor corrections to the evaluation of the time 
scales of the phenomena, but — because at the beginning 
of the simulations the chains are locally straight, and be- 
cause after the collapse the hydrodynamic interactions 
would be dominated by the attractive forces — it would 
not change the kinetic pathways and the relative stability 
of the intermediate states. 

3. RESULTS 

This section is organized as follows: Section lTTI reports 
qualitatively the results of our simulations, describing the 
various possible intermediate configurations and the typ- 
ical pathways. Section presents a more quantitative 
description of the collapsed states in term of the corre- 
lation function between tangent vectors along the chain 
and of the energy levels associated with each different 
collapsed state, which confirm that the toroid is the sta- 
ble configuration, and show how the collapse proceeds 
through a cascade of subsequently more energetic favor- 
able metastable states — the various multiple headed 
racquets shapes, as predicted from theory |23|. The de- 
cay rate of uncollapsed states is investigated as a func- 
tion of a proper combination of L, Lq and Lp in Section 
13.31 together with a statistical analysis of the preferred 
collapsed configurations. Section provides a compari- 
son between the results obtained with the different com- 
putational models (phantom chain and bead- rod chain). 
Some preliminary results of the effect of shear fiow on 
the kinetics and pathways of collapse are presented in 
Section [3. 51 



3.1. Molecule configurations and collapse patliways 

Figure Q] shows a typical snapshot of an ensemble of 
semiflexible chains during collapse. The 3-D structures 
of the molecules are shown here as 2-D projections on 
a convenient plane. At a given time molecules coexist 
in different configurations, from completely uncollapsed 
(molecule (a)), to partially collapsed (molecules (b), (f) 
and (g)), to higher order racquet (molecules (c) and (e)), 
to toroid (molecule (d)). 

The fraction of collapsed molecules increases monoton- 
ically with time, confirming that collapsed states are en- 
ergetically favorable in a bad solvent. Many other en- 
sembles with different values of L/Lp and Lg/Lp show a 
similar behavior. However, if the attractive forces are not 
strong enough or the bending stiffness is too high there 
is no evidence of collapse and the open conformations 
persist in time. As already observed in [Tsf . a sequence 
of collapse can evolve following two different generic pat- 
terns. The first pattern involves the direct formation of a 
torus from the open conformation, when the two ends of 
the chain meet while the tangent vectors are locally par- 
allel. In this instance, the chain collapses without forma- 
tion of any intermediate, the end to end distance drops 
quickly to a value of the order of the torus dimension, 
and then fluctuates as the chain keeps folding into the 
final structure. Meanwhile the total energy of the chain 
keeps decreasing monotonically till reaching the equilib- 
rium value. Conversely, the second pattern involves the 
formation of intermediate metastable states: initially a 
single racquet-headed shape which folds quickly into a 
short-lived, 3-heads racquet configuration, and then re- 
arranges as a toroid via a subsequent folding. Figure El 
shows the two different pathways. We notice here that 
the direct formation of a torus is a more unlikely pattern 
in 3-D — and in fact observed with scarce frequency — 
than in 2-D. Whereas in 2-D the relative orientation of 
the two ends is described by a single angle, two angles 
are required in 3-D and both need to be close to 2% for 
the direct formation of a torus to take place. Most of the 
chain will therefore initially form a single headed racquet, 
as described below ^ . 

The time evolution of the total energy of an ensemble of 
collapsing chains is shown in Figure 13 The total energy 
U of the chain during collapse is defined as the sum of the 
bending energy f/bend and the interaction energy along 
the chain UsoW- In the discretized form, 

W-2 N N 

U ^ kTLp cos^ e, + Y, "=°iv,ifc' (7) 

i—l j—1 k—j + 1 

where 9i is the bond angle and Wsoivjfe is the interaction 
energy evaluated between the beads j and k along the 
chain. 

The time series in Figure|21reports most of the possible 
evolutions observed in our simulations. Few molecules re- 



4 



formation of torus through 
higher order racquet 




direct formation of torus 



100 



150 



200 



250 



20 
& PM5 
_b ^ 10 

° J. 5 




"slow" raveling of the chain 



sudden drops — chain folding on itself 



50 100 150 

time (unit of t^tQ 



200 



250 



FIG. 2: Comparison between time series of end to end dis- 
tance K and total energy for direct formation of toroid and 
formation through intermediate states. 



main in the uncoUapsed states for the whole duration of 
the run, equivalent to 250 Tdiff^L . For these molecules, 
the thermal fluctuations never lead two branches of the 
chain close enough together to initiate the collapse. The 
other molecules typically reach lower energy states by col- 
lapsing in sequentially higher order metastable conforma- 
tions. The energy levels for each racquet-head shape of 
increasing order are well separated, and there is a rather 
large gap from the highest order racquet observed (6* ^ 
order) and the stable tori, as expected from theory |23l| . 
Most of the chains initially form a single head racquet: 
after the first monomer- monomer contact and the forma- 
tion of the racquet head, the two sides of the chain tend to 
become parallel and form a neck region, which ensure the 
maximum number of monomer- monomer contacts with- 
out increasing the bending penalty. The molecules then 
fold again on themselves, forming higher order shapes. 
Depending on the way in which the molecule folds back 
on itself, the next configuration can be a 2-heads or a 
3-heads racquet shape. Figure 13 shows examples of both 
pathways. Successive folding leads to higher order rac- 
quets. The transition between each single metastable 
conformation is very rapid and driven by the determin- 
istic attractive and bending forces. In contrast, the initi- 
ating event for each successive collapse is related to the 
Brownian motion of the molecule, and therefore the time 
interval between folding events is extremely variable, and 
completely random for the simulations that we have per- 
formed. In the simulations of longer chains ( L = 10 Lp ) 
we observe combinations of the conformational elements 
described above, e.g., chains partially folded in a racquet 
shape at one end while still uncoUapsed at the other. 
Those conformations have intermediate energies between 
the well-defined levels corresponding to racquets and tori 
and quickly disappear in favor of more compact and or- 
dered structures. In contrast, the metastable high order 



racquets are extremely long-lived, and the energetic bar- 
rier for reaching the equilibrium shape is of the order of 
several k^T . While the partial unfolding of a torus in a 
racquet is never observed, simulations show the forma- 
tion of a torus via metastable configurations, and this, 
together with the lower total energy shown by the torus 
(see Figure O), confirms once more that the latter is the 
equilibrium configuration. 



3.2. Quantitative description of collapsed states 

The configurations and the shapes of semiflexible 
molecules in bad solvent can be studied through direct vi- 
sualization of the collapsing chain. While this method is 
a valuable way to gather information on the kinetics and 
pathways of collapse, it is time consuming and not quan- 
titative. We therefore define a spatial correlation matrix 
M to analyze in a more systematic way the collapsed 
shapes. The elements of M are defined as My = • Uj. 
With this definition, M is symmetric; a perfectly straight 
rod has A/y = 1 Vz,j, as all the connecting normalized 
vectors have the same direction and orientation. An un- 
coUapsed semifiexible chain with persistence length Lp 
will have Ma = 1, while the off-diagonal terms My would 
be smaller and decaying exponentially with the distance 
along the chain between the connectors i and j. 

When the chain collapses, the actual shape assumed by 
the molecule can be inferred from the spatial correlation 
between connectors. The spatial correlation matrix for a 
single head racquet, for example, will be a block matrix: 
one block will show values close to 1, corresponding to 
the almost straight filament in the neck of the racquet, 
followed by a region with rapidly decreasing values from 
1 to —1 — the racquet head — and by another block 
with values close to —1, indicating the other filament of 
the neck , with same direction but opposite orientation. 

The spatial correlation matrix can be simplified into a 
spatial covariogram for the molecule. For standard sta- 
tistical analysis, the covariogram is defined as 

N N 

c{h) = l/N{h)Y,T.i^^-^'}i^J-^'} (8) 



where h is the distance (in time or in space) between 
the observations Zi and Zj, N{h) is the number of ob- 
servations at a distance h and fi is the expected average 
value of the observation. In the present analysis, for each 
molecule we construct the spatial covariogram Csp assum- 
ing Zi = Ui and /i = 0. We therefore obtain a scalar func- 
tion of the distance along the chain s, bounded between 
+1 and —1, which contains all the relevant information 
about the shape of the molecule. 

Figure 21 shows three relevant examples of this type of 
analysis: the actual shape of the semiflexible molecules 
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FIG. 3: Time series of total energy for an ensemble of molecules, with labels of the corresponding configurations and actual 
shape of (a) 1 head (b) 2 heads (c) 5 heads racquets and (d) torus. 



is shown in comparison with the spatial correlation ma- 
trix and the spatial covariogram. For a uncoUapsed chain 
(Figure ^a), the matrix shows randomly alternating re- 
gions of positive and negative correlation along the length 
of the chain, corresponding to the typical open configu- 
ration of a semiflexible chain in good solvent. The cor- 
responding covariogram shows how the connectors cor- 
relation decays from 1 to values around along the 
length of the chain; the decay is exponential at short 
distances and related to the stiffness of the chain, i.e., 
Csp — exp(— s/Lp). At longer distances, the correlation 
randomly oscillates between negative and positive values. 

For a racquet shape with multiple heads, M has a 
structure of alternating blocks with strong positive and 
strong negative correlation, as shown in Figure ^J-b. The 
5 straight filaments in the neck of the racquet correspond 
to these 5 blocks in the matrix; the size of the blocks is 
approximately the same and indicates the length of the 
racquet neck. The covariogram shows 4 zero-crossing, 
each corresponding to a head of the racquet, and can be 
fitted quite accurately with a 4*^ order harmonic; the pa- 
rameters of this fit are related to the curvature of the 
racquet heads. 

For a torus, the spatial correlation matrix appears 
completely different (Figure EJc), showing a structure 
with diagonal bands. In fact, the correlation between 
pairs of connectors at the same distance is the same, re- 
gardless of their position along the chain; this translates 
into bands along the diagonals of the matrix. The dis- 
tance between a diagonal with values close to 1 and one 
with values close to — 1 corresponds to the diameter of the 



collapsed torus. The same information can be extracted 
from the covariogram, which is a perfect sinusoid: the 
period of the sinusoidal function corresponds to tt times 
the diameter of the torus, and the number of complete 
periods in along the chain indicates the number of loops 
formed by the collapsed molecule. 



3.3. Kinetics of collapse 

The fraction of collapsed semiflexible molecules in a 
bad solvent increases monotonically with time, as already 
stated. We have studied the effect of the relevant length 
scales of the phenomenon, Lp, Lq and L, on the kinet- 
ics of this collapse. By systematically varying the to- 
tal length of the chain (L = 3Lp, 5Lp, 8Lp and lOLp), 
as well as the strength (uq — Q.2bkT/a?, Q.bkT/a? and 
\kT/a?) and the range (i?attr = 5/AOLp and l/6Lp) of 
the attractive forces, we have explored the effect of these 
parameters on the decay rate of uncoUapsed molecules. 
We define the time to collapse tcoii as the time required by 
the molecule to form the first metastable collapsed struc- 
ture. We obtain from simulations the average value of tcoii 
in an ensemble of collapsing chain under each condition, 
normalized with Td\ff,Lj,- The inverse of icoii is taken as a 
measure of the decay rate Dcoii of uncoUapsed molecules. 
The distribution of the times to collapse follows the ex- 
pected lognormal distribution, typical of such dynamical 
processes, with a peak close to the average value and a 
long tail. 

At a given value of L/ Lp, the decay rate increases. 
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FIG. 4; (a) Actual shape, (b) correlation matrix and (c) spatial covariogram (with fitting curve) for (1) an uncoUapsed molecule, 
(2) a metastable racquet-head shape and (3) a stable torus. 
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FIG. 5: Masterplot of the decay rate of uncoUapsed molecules. 



i.e., the molecules collapse faster, for decreasing values 
of Lqj i-G., for increasing values of uq. Similarly, the 
decay rate increases for longer chains at a given value of 
the attractive potential and persistence length. We can 
collapse the values of Dcoii on a single curve as a function 
of a proper combination of the two dimensionless ratios, 



L/Lp and L^/Lp, {L/Lp)^^^ {Lp/Lo). Figure[Slshows the 
mastercurve for the decay rate. 



Table reports the statistics of chains conformations 
for different values of L/Lp at two different times, i.e., 
tl = 100rdiff,Lp and t2 = 200Td\ff,L^- At tl, the the single 
head racquet is a conformation much more likely than the 
torus for all the explored values of L/Lp. This confirms 
that the direct formation of a torus is an unlikely event, 
as discussed in Section [3.11 and most of the chains col- 
lapse through a sequence of folding events. At t2, the per- 
centage of molecules in the 1-head racquet configuration 
does not vary substantially with respect to tl, because 
— while uncoUapsed molecules fold to form new 1-head 
racquets — some of these configurations further collapse 
in higher order metastable racquets and stable tori. The 
number of chains collapsed in tori or multiple head rac- 
quets monotonically increases within the observed span 
of time. For shorter chains the likelihood of tori and mul- 
tiple racquets is similar, while the metastable states seem 
to be initially more favorable for longer chains {L/Lp = 8 
and L/Lp = 10). Overall, the results show that the inter- 
mediates conformations in the collapse phenomenon are 
extremely long-lived: only 8.4% of the simulated chain 
has reached the equilibrium, torus state after 200rdiff.L . 
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Conformation 


L/Lp = 3 


L/Lp = 5 


L/Lp = 8 


L/Lp = 10 


tl 


12 


11 


12 


11 


12 


11 


12 


open chain 


62.0 


43.3 


48.0 


30.0 


43.3 


36.7 


49.0 


36.8 


single-head 


31.0 


43.3 


36.0 


42.5 


33.4 


30.0 


24.5 


23.7 


muhiple-heads 


3.5 


6.7 


8.0 


15.0 


20.0 


26.7 


20.4 


31.6 


torus 


3.5 


7.9 


8.0 


12.5 


3.4 


6.7 


6.1 


7.9 



TABLE I: Statislics ol chain conformations for different L/Lp 
at 11 = 100rdiff,Lp and 12 = 200rdiff,Lp (% values). 

3.4. Comparison of different computational models 

We have performed simulation of collapsing chains 
with different computational models, in order to verify 
the results obtained and the appropriateness of the model 
chosen. In particular, we have performed simulations 
with a bead-rod model, where the chain connectors are 
rigid rods rather than stiff elastic spring. As noted in Sec- 
tional the bead-rod model is a constrained model which 
ensures local inextensibility of the chain, but is computa- 
tionally more expensive, due to more stringent time step 
requirements for convergence. While local inextensibility 
is fundamental when analyzing the polymer contribution 
to the stress tensor its effect on the pathways and 
dynamics of collapse is less relevant, as confirmed by our 
results. In fact, we observe qualitatively similar path- 
ways to collapse in the simulations performed with the 
bead-rod model. The preferred pathway includes the for- 
mation of metastable racquet shapes, similarly to what 
observed in the bead-spring simulations. The kinetics of 
the collapse are slower for the locally inextensible model; 
we believe that this is due to how the constraints are 
imposed, i.e., by projecting the random forces onto the 
constraints, and therefore not allowing any fluctuations 
along the chain. A more detailed comparison between 
the collapse kinetics two model is prevented by the high 
computational cost of simulating the constrained system. 

We also perform simulations with a phantom bead- 
spring chain, i.e., a chain with no repulsive hard core. 
Once more, we find that the dynamics of collapse are 
qualitatively the same, although the actual shapes of the 
torii and of the racquets formed by the phantom chain 
are slightly different. The absence of hard core repulsion 
allows in fact for tighter structures, as shown in Figure 
El In particular, in the metastable multiple racquets the 
neck for phantom chains appears much narrower than for 
the chain where excluded volume interactions are taken 
into account, while it is noticeable that the shape and the 
curvature of the racquet head appears unchanged. The 
traces in time of the energy level for the phantom chain 
show the features already discussed in Section mi well 
defined energy states corresponding to metastable struc- 
tures, with rapid transitions between them (see Figure 
El inset). However, the average life of metastable states 




FIG. 6: (a) Metastable racquet-headed stale for a phantom 
chain and (b) corresponding correlation matrix. Inset in (a): 
energy traces in time for collapsing phantom chains. 

appear slightly shorter than in the complete model: in- 
tuitively, this is due to the wider range of motion of the 
phantom chain, which permits folding paths otherwise 
prevented by the excluded volume interactions. 



3.5. Effect of shear flow 

We discuss here the results of some preliminary work 
on the effects of steady shear flow on the collapse dynam- 
ics of semiflexible chains. In particular we monitor the 
decay rate Dcoii at different flow strengths, and we com- 
pare the pathways of collapse with those at equilibrium. 
Flexible chains expand in shear flow, i.e., the average 
end to end distance increases while the molecules tend 
to orient with the flow; the effect of the flow is there- 
fore counteracting that of the attractive forces, which 
tend to form compact globula. Stiff semiflexible chains, 
with Lp > L shrink in shear flow, due to a buckling 
instability j25||; such rigid chain however would not col- 
lapse to form tori and multiple racquets in bad solvent, 
as the open, rod-like shape is stable at equilibrium [23j . 
In this work, we consider the collapse dynamics of less 
stiff semiflexible molecules: is not clear a priori what to 
expect for the collapse dynamics under shear flow. We 
show here that strong enough shear flow speeds up sen- 
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FIG. 7: Plot of the decay rate Dcoii vs. Wi for L/Lp = 3, 5, 
8, uo = IfcT and R^tt, = 10/3. 




FIG. 8: Collapsed configurations of semiflexible chains in 
shear flow at Wi = 10; (a) L/Lp = 8, (b) L/Lp = 5. 



sibly the collapse kinetics, increasing the likelihood for 
the two ends of the chain to meet. Figure [7| show the 
trend of the decay rate as a function of Wi, where Wi 
= 7Tdiff,Lj,, for 3 different values of L/Lp. At Wi = 0.1, 
the decay rate does not change much with respect to the 
equilibrium value, and there is not a clear trend among 
the different L/Lp. At Wi = 1, the increase in decay 
rate becomes instead significant, as the average time to 
first contact halves; for Wi = 10, the effect of shear flow 
is completely dominating and the decay rate under the 
3 different conditions is almost equivalent. Under shear 
flow, the semiflexible molecules undergo a sequence of 
compressions and extensions while tumbling in the plane 
of shear; in the compression region, the molecule tends to 
shrink, therefore raising the probability that two sections 
far apart along the chain backbone come in contact. The 
attractive forces then trap the molecule in the partially 
collapsed configuration, preventing it from straightening 
back in the extension region. However, this continuous 



tumbling of the chain and alternation of extension and 
compression does not allow the formation of the higher 
order, compact structures observed at equilibrium. In 
particular, the shear flow seems to inhibit the formation 
of the torus, which was never observed in the high shear 
flow simulations. Figure |S| shows two typical configura- 
tions observed in strong shear flow; while they show the 
basic features of the metastable racquets, the structures 
are less compact and do not reach the same level of en- 
ergy of the equilibrium conditions. 

4. DISCUSSION 

The dynamics of collapse of semiflexible molecules in 
bad solvent has been investigated in the literature. Di- 
rect, visual experimental evidence of the formation of 
torii and high order racquets has been obtained study- 
ing biop olymers with sufficiently large persistence lengths 
13, 1141 . Theoretical work and Monte Carlo simulations 
llL I12L l23j have confirmed that the torus is the sta- 
ble conformation, but intermediate, metastable racquet 
shaped conformations exist and are long lived. Previous 
Brownian Dynamics simulations in 2-D have firstly 
suggested the possible general pathways to the dynamics 
of collapse. 

The present work provides new insights on the collapse 
of semiflexible chains, confirming the suggested pathways 
through 3-D simulations, and showing the effect of the 
two relevant dimensionless ratios L/Lp and L^/Lp on 
the kinetics of collapse. The time series of the total en- 
ergy of each chain show clearly the fast transition be- 
tween well defined energy levels, corresponding to differ- 
ent metastable intermediate shapes, and the final, sta- 
ble torus. We also showed that two possible pathways 
towards the stable conformation are possible: direct for- 
mation of a torus and successive folding of the chain in 
progressively higher order racquets. The latter one ap- 
pears to be more probable for the range of parameters 
investigated. The decay rate of uncoUapsed states for 
all the different conditions of the simulation can be plot- 
ted as a function of a proper combination of L/Lp and 
L^/Lp. Further work is needed to confirm the validity 
of this empirical scaling and verify it in a wider range of 
parameters. 

We also performed simulations with different compu- 
tational models, therefore confirming that the observed 
metastable conformations and pathways are general and 
do not depend on the choice of the model. However, the 
kinetics, i.e., the actual value of the decay rate can be 
different: in particular, for a constrained bead-rod model 
the kinetics of collapse appear slower, possibly due to the 
way the Brownian noise is projected onto the constraints. 

We show preliminary results on the kinetics and path- 
ways of collapse of semiflexible chains in shear flow. 
Strong enough shear flow increases the decay rate of un- 
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collapsed states, but prevents the formation of the com- 
pact, well defined configuration observed at rest, and 
seems to inhibit the formation of the torus. While these 
results have been obtained in unbounded shear flow, fu- 
ture work will consider the interaction of the collapsing 
molecules with a solid surface, which could be attractive 
or repulsive for the chain. Also, future possible direction 
of research include the effect of local defects along the 
chain, which make the molecule locally more flexible or 
stiffer. 
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